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We analyse the ground state of a particle in a double-well potential, with a cylindrical symmetry, 
and in the presence of a magnetic field. We find that the azimuthal quantum number m takes the 
values m = 0, 1, 2 . . . when we increase the magnetic field. At critical values of the magnetic field, 
the ground state is twice degenerate. The magnetisation shows an oscillatory behaviour and jumps 
at critical values of the magnetic field. This phase transition could be seen in the condensate of a 
dilute gas of charged bosons confined by a double-well potential. 



I. INTRODUCTION 

It is a well-known factlll that in the absence of a mag- 
netic field the ground state of bosons is non degener- 
ate, and therefore has the symmetry of the hamiltonian. 
Mathematically this result from the fact that the kernel 
of the operator e~*^ is positive. This last property no 
more holds in the presence of a magnetic field so that de- 
generacy of the ground state may be expected, as well as 
symmetry breaking in it. One-body systems already mas4 
show this phenomenon. Indeed Lavine and O'Carrolcl 
proved the existence of spherically symmetric potentials 
for which, in the presence of a magnetic field, the ground 
state has a non-vanishing value for the z component of 
angular momentum, so that the rotational symmetry is 
broken. 

Further examples were provided by Avron, Herbst and 
SimonBu. On the opposite side, these last authors were 
able to prove that for the hydrogen atom the symmetry 
is not broken, as well as in the case where the poten- 
tial is monotonically increasing with the distance. These 
authors, however, mainly concerned with problems of 
atomic physics, did not discuss the degeneracy and the 
physical significance of it. 

On the other hand, two of us analysing the problem of 
a particle confined to a disc or an annulus in the pres- 
ence of a magnetic field found that the ground state was 
degenerate in the case of an annulus and for a disc with 
Neumann boundary conditions (with Dirichtlet boundary 
conditions in the disc case the degeneracy disappears.) 
The degeneracy appears each time the magnetic field 
reaches a critical value and the magnetisatioiinjumps at 
these critical values, which form a discrete setu. 

Motivated by these results we consider in this arti- 
cle a class of systems for which similar phenomena oc- 
cur. Namely we analyse the ground state of a particle 
in three dimensions moving in a double-well type poten- 
tial, cylindrically symmetric, and submitted to a constant 
magnetic field in the z direction. 



We find that the ground state has an azimuthal mo- 
mentum hm taking increasing values m — 0,1,2, ... when 
we increase the magnetic field B. At critical values of B 
(Bm) the ground state is twice degenerate between the m 
and the m + 1 state. Moreover the magnetisation jumps 
at these critical values and shows in general an oscilla- 
tory behaviour reminiscent of the well known de Haas 
von Halphen oscillations in solid state physics. 

We show that this phenomenon can be understood by 
an analysis of the minima of the potential energy, fixing 
however the angular momentum to its quantised value 
hm. In the two-dimensional case we can use the WKB 
method and obtain bounds on the energy in order to 
estimate the critical fields. But in general, we had to 
compute numerically the energies and compare them to 
estimates based on trial wave functions. The agreement 
is quite good in general. 

Concerning possible experimental verifications of these 
effects, which require basically to have a potential which 
has a minimum sufficiently far from the origin, we could 
think of two cases. The first one would be in some 
molecules where proton dynamics could be described 
by such an effective potential. The second one, more 
thrilling, would be the case of charged bosons undertak- 
ing a Bose-Einstein condensation. Our results suggest 
that in this case, the bosons would undertake a phase 
transition in their condensate, when we apply an in- 
creasing magnetic field. This phase transition would 
manifest itself by appearance of oscillations in the mag- 
netisation, which would jump at certain critical values 
of the magnetic field. 

II. THE MODEL 

We will consider the case of a particle of mass fj,, charge 
q, in a potential V with a cylindrical symmetry, submit- 
ted to a magnetic field B in the z direction. We do not 
consider the effect of the spin of the particle. We choose 



1 



for a unit of energy Vq, and length vq, both being charac- 
teristic of the potential. The dimensionless hamiltonian 
reads if r = 

(ieV-l)2 + y(r,z) (1) 

where 

' (2) 



measures the importance of the quantum effects and the 
vector potential in the symmetric gauge is given by 
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B — — " B being the dimensionless magnetic field. 
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Thanks to the cylindrical symmetry, we can replace 
the z component of the angular momentum Lz by its 
eigenvalue em so that the reduced hamiltonian reads 
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The ground state energy of this hamiltonian and the cor- 
responding eigenfunction will be denoted Em and ipm- 

It remains to specify V. We will basically consider a 
double-well potential of the form: 



V{r, z) 



r'* + z^ 



,2 Ji 



2(r^ + z^) +vr^z 
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with V satisfying v > —2, so that V is bounded from 
below. If V is equal to we can decouple the motion in 
the z direction form the one in the plane perpendicular 
to the magnetic field. This is what we will call the two- 
dimensional case. If w = 2, we have in three dimensions 
a potential with spherical symmetry. 

We have chosen this double- well form because if we had 
taken the simple well V ^r"^ + z"^ + 2(r^ -I- z^) -I- vr'^z^ 
with w > it follows fromO that the ground state is not 
degenerate and corresponds to to = 0. 

A physical quantity of interest is the magnetisation 
in the ground state 



M = 



dE 
"dB 



(6) 



in units ^ro^^ 

We will denote by the ground state energy of the 
hamiltonian 



with 
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+ V„,{r,z) (7) 
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and by 

E„i = e„ - emB (9) 

the ground state energy of i?™ given in so that the 
real ground state energy is given by 

E = inf E,n (10) 

since obviously negative to give a larger energy. 

Finally we will use the following useful scaling property 
of the energy e„i 



where 



V s > 



(11) 



(12) 



is the parameter multiplying in the potential. Equa- 
tion {nW follows simply from the scaling transformation 
: sr'^ and — ^ sz^. This relation shows that we 

have effectively a two parameter dependence of the en- 
ergy Cm in general and a one parameter dependence in 
the two dimensional case. 

The choice s = |A| or s = tos (to > 1) shows that 
large magnetic field or large angular momenta correspond 
to the semi-classical limit. In fact we shall see that in 
the classical limit e — > ground state with to 7^ are 
favoured inducing ground state degeneracies at some val- 
ues of the magnetic field. It thus appears that the ten- 
dency to have a ground state with the same symmetry as 
the hamiltonian and therefore non degenerate is an effect 
due to quantum mechanics. 



III. THE CLASSICAL LIMIT 

One can gain some qualitative understanding of the 
problem by looking at the classical limit of it. This means 
that we neglect the quantum kinetic energy and define 
the ground state energy as 



E = inf inf \ym{r, z) — emB] 

m>0 (r.z) 



(13) 



where 
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' 7^ + — r^ -~2r^ + v r^z"^ (14) 



and consider that to is an integer. 

Two cases need to be considered separately: < 2 
and V >2. li\v\ < 2 we denote by Xm and tm respectively 
the values of and z^ which minimise the potential V^m 
and we find 



^rn. — 1 



(15) 
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On the other hand, considering for a while m as a con- 
tinuous variable, the absolute minimum of Vm — emB is 
given by 



em 



B 



(16) 



From (|15|) this gives an absolute minimum of Vm — emB 
given by 



Xrh 



m — 



1 



and therefore 
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(17) 
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In considering the variable m as a continuous one we 
have treated the problem purely classically and the cor- 
responding "ground state" energy is 
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We know that m is a discrete variable but for consistency 
we must consider e as a small number. Then if m des- 
ignates the integer part of m, we have m = m + and 
if < < ^, the ground state has the quantum number 
TO, whereas if ^ < < 1 it has to -f 1. 

From this analysis we conclude that if Bm-i < B < 
Bm where 



B„, = e(l + -)(2TO 



1) 



(20) 



the ground state has the quantum number m. Hence we 
see that by increasing the magnetic field, we find in in- 
creasing order the values of to = 0, 1, 2, ... and an infinite 
set of critical values of the magnetic field exist, Bm 
for which the ground state is twice degenerate, being 
both TO and to + 1 . 

This picture is entirely confirmed by the numerical re- 
sults in the quantum case. It is also quite interesting to 
look at the magnetisation. In the state whose quantum 
number is to, we have 



M„ 



so that using (|T^) 
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when Bm-i < B < Bm- 

This shows that the magnetisation has an " oscillatory" 
type of behaviour reminiscent of the familiar de Haas von 
Halphen one in solid state physics and that the magneti- 
sation jumps at the critical values of the magnetic 
field, the jump being given by 



AM„ 



1-! 



(23) 



Once again this general behaviour is reproduced by 
the numerical results in the quantum case and the spac- 
ing between the values of the critical field is rather well 
represented by formula (|20| ) when to > 1. In the two- 
dimensional case, i.e. z; = and neglecting the trivial 
z dependence, we can proceed further and look at a re- 
ally semi-classical approximation namely WKB, for the 
ground state energy 



/ dr^e 



Vm{r) 



where 
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and the ground state energy is Em = &m — emB. 

In fact this WKB approximation will give the best an- 
alytical results, apart from the variational estimates for 
the energy, which give unfortunately only exact upper 
bounds on the energy. 

When the potential has spherical symmetry v = 2, 
quantum effects are much more important and the clas- 
sical analysis gives only that the ground state has to = 
if _B < 2e, is degenerate between to = and to = 1 when 
2e < B < 4e, has possibly to = 0, 1, 2 for 4e < B < 6e 
and so on. This only suggests that we have again the 
increasing sequence of to, when we increase the magnetic 
field and that critical values appear near 2eTO. 

When w > 2, we find that to = is the ground 
state except when B = 2eTO, where it is degenerate be- 
tween TO and 0. We may note however that the classical 
ground state correspond to points (r — 0,z — ±1) in 
configuration space for to = 0, whereas it corresponds 
to two circles (r 



2B ' 



2B 



for TO = 1 and 

2e < B < 4e, so that the wave function can be more 
spread in the m = 1 state than the in the m — state, 
and that the kinetic energy of the m — I state is lower, 
favouring the to = 1 state. Hence we should expect, at 
least when e is small, a ground state with m = for small 
fields and a ground state with to = 1, when 2e < B < 4e. 
A similar argument can be given for the higher values of 

TO. 

Finally, it is worth noticing that if we had taken a 
simple well type potential 



Vir,z) 



r'^ + z^ + 2{r^ + z^) + V z^ 



(26) 



the classical analysis gives a ground state with m = 0, at 
least when w > — 1. This is a correct result when v > 
at the quantum level. 



3 



IV. NUMERICAL RESULTS AND VARIATIONAL 
BOUNDS 



It is quite useful to undertake a numerical analysis of 
this problem. We have used a finite element method, 
choosing for the basis a product of two triangles func- 
tions. We discuss separately the two-dimensional prob- 
lem and the three dimensional ones. 



A. Two dimensions 






simulation 


* 


estimation 



We first give pictures of the ground state energy for 
two typical values of e, a small (e = 0.03) and a large one 
(e — 0.5) as a function of the magnetic field B. (figure 
The cusps at the critical values of B indicate a jump 
of the corresponding magnetisation. 
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FIG. 1. Energie and magnetisation dependence of B for 
e = 0.03 and e = 0.5 



This last quantity shows first a diamagnetic behaviour 
at small field, but then a paramagnetic - diamagnetic 
oscillation at least when e < 0.3. Beyond this value the 
magnetisation is entirely negative (figure |l| bottom right). 
We can also note that when B becomes large the mag- 
netisation tends to — e, its value in the Landau regime. 

The results clearly indicate that we go progressively 
from the states with m = 0,1,2... by increasing the 
magnetic field and that the magnetisation jumps at the 
critical values. The effect is more pronounced in the clas- 
sical regime. All these results are in qualitative argument 
with the classical picture presented before and the agree- 
ment is even quantitative when e — 0.03 for example. 



FIG. 2. Comparison of the jumps given by simulation and 
the estimation 



The jumps of the magnetisation given by formula (|22^^ 
are reproduced (figure 0) with a precision of less than 
1 percent when e — 0.03, and the spacing between the 
critical values of the magnetic field 



Bri 



Bn 



= 2 + A„ 



(27) 



is given by A^, < 0.04 if m > 1 and e = 0.1. A^, de- 
creases when m increases in agreement with the scaling 
relation Bm = {2m + l)e, so that the simple classical for- 
mula reproduces rather well the results. By contrast, the 
jump between the m = and the rn ~ 1 state is largely of 
quantum mechanical origin, as well as the precise values 
of the critical fields. 




FIG. 3. AI values of the ground state depending on B and e 



Figure || describes the various regions in the e— B plan. 
We can note that even when e > 0.25 a linear relation 
exists between Bm and e, as in the classical regime, which 
is a bit surprising. 
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m = B = 0.34 E = 0.2 m = 1 B = 0.34 E = 0.2 
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Variat . 
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i-z 
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U.DO 


-U.y4Uo 


n 07 
V.Zl 


2-3 


-0.9403 




0.64 


-0.9401 


0.28 


3-4 


-0.9401 


-0.9397 


0.63 


-0.9399 


0.29 


4-5 


-0.9399 


-0.9395 


0.61 


-0.9397 


0.30 


5-6 


-0.9396 


-0.9392 


0.59 


-0.9394 


0.30 


6-7 


-0.9392 


-0.9389 


0.56 


-0.9390 


0.30 



TABLE II. Energies Em at the seven first degeneracies 
with e = 0.03 




FIG. 4. Potential and eigenfunction of the ground state for 
m — (left) and m — 1 (right) with Vm ( — ) ~ orbital kinetic 
energy (— • — ) + double- well ( ) and energy levels (• • ■) 



It is also interesting to look at the eigenfunctions when 
the magnetic field reaches its critical value. In figure ^ 
we give pictures of them at the critical value between the 
state m = and m = 1 when e = 0.2. We see that 
their maxima are located very near the minimum of the 
potential. 

Finally we compare the results with two theoretical 
estimates: first of all the WKB one, and a variational 
one. This last estimate is based on the following two 
parameters trial wave function 



(28) 



The variational upper bound on the energy can be ex- 
pressed in terms of Weber cylindrical functions, but we 
directly computed the corresponding integrals. 



Deg. 


Simul. 


WKB 


<5 % 


Variat. 


5 % 


0-1 


0.0313 


0.0314 


0.23 


0.0317 


1.16 


1-2 


0.0944 


0.0942 


-0.15 


0.0946 


0.23 


2-3 


0.1573 


0.1571 


-0.13 


0.1574 


0.09 


3-4 


0.2201 


0.2198 


-0.13 


0.2203 


0.06 


4-5 


0.2830 


0.2826 


-0.12 


0.2830 


0.02 


5-6 


0.3457 


0.3453 


-0.12 


0.3458 


0.02 


6-7 


0.4085 


0.4080 


-0.12 


0.4085 


0.00 



TABLE I. Magnetic field _B„ 
with e = 0.03 



at the seven first degeneracies 



Tables |,|I,[II and IV give a comparison of the results 
for two values of the parameter e, and for the critical 
fields. Excellent agreement is found for the variational 
method (maximal error of the order of 2 % when e = 0.5). 
WKB works quite well when e is small (e = 0.03) as 
expected, but even better on the energies when e = 0.5 
and the error does not exceed 1%. 



Deg. 


Simul. 


WKB 


6 % 


Variat. 6 % 


0-1 


1.538 


1.661 


7.95 


1.508 


-1.98 


1-2 


2.747 


2.811 


2.33 


2.743 


-0.15 


2-3 


3.842 


3.882 


1.06 


3.842 


0.02 


3-4 


4.891 


4.919 


0.56 


4.894 


0.05 


4-5 


5.920 


5.940 


0.34 


5.924 


0.07 


5-6 


6.941 


6.954 


0.18 


6.943 


0.02 


6-7 


7.953 


7.964 


0.12 


7.956 


0.02 



TABLE III. Magnetic field Bm at the seven first degenera- 
cies with e = 0.5 



Deg. 


Simul. 


WKB 


8% 


Variat. 5 % 


0-1 


0.220 


0.232 


0.97 




0.55 


1-2 


0.685 


0.686 


0.04 


0.690 


0.25 


2-3 


1.159 


1.159 


-0.02 


1.163 


0.16 


3-4 


1.639 


1.638 


-0.03 


1.642 


0.12 


4-5 


2.122 


2.122 


0.00 


2.125 


0.12 


5-6 


2.609 


2.608 


-0.02 


2.612 


0.07 


6-7 


3.098 


3.098 


0.00 


3.101 


0.07 



TABLE IV. 
with e = 0.5 



Energies Em at the seven first degeneracies 
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B. Three dimensions 



For the spherically symmetric potential (t; = 2), figure 
^ gives the ground energies a well as the corresponding 
magnetisation for two different values of e : 0.03, 0.5. 





FIG. 5. Energie and magnetisation dependence of B for 
e = 0.03 and e = 0.5 



between the critical values of B predicted by the crude 
classical estimate: 



(29) 



is satisfied with a precision of 25% at m = 1 and becomes 
more accurate when m increases, at least in the range 
e < 0.1. 

Our best variational estimate for the energy was made 
with a three parameter trial wave function 



V'a./s.c = r'"e-"'' -/3(vV^-C) (30) 

Table gives the values of the critical field Bm and Table 
VI the corresponding ground state energies, when e = 
0.05 estimated by the variational method and computed 
with the simulation. 





Simulation 




Variational 




Dcg. 


Bni 


Em 


B,n 


Sb % 


Em 


Se % 


0-1 


0.1180 


0.1206 


2.17 


-0.8986 


-0.8982 


0.38 


1-2 


0.2381 


0.2310 


-2.94 


-0.8966 


-0.8966 


-0.01 


2-3 


0.3549 


0.3509 


-1.15 


-0.8946 


-0.8946 


-0.06 


3-4 


0.4686 


0.4616 


-1.49 


-0.8925 


-0.8925 


-0.00 


4-5 


0.5829 


0.5785 


-0.74 


-0.8901 


-0.8901 


-0.01 


5-6 


0.6961 


0.6905 


-0.80 


-0.8876 


-0.8876 


-0.00 



TABLE V. Magnetic field Bm and energies Em at the six 
first degeneracies at e = 0.05 



Once again we see that the values of m in the ground 
state increases with B, and that the magnetisation jumps 
at critical values Bm of the magnetic field, where the 
ground state is doubly degenerate. These results are in 
qualitative agreement with the classical analysis. Figure 
y summaries the results in the e - B plane. Notice that 
in this 




B 

FIG. 6. M values of the ground state depending on B and e 





Simulation 




Variational 




Deg. 


Bm 




Bm 


Sb % 




Se % 


0-1 


2.7576 


0.9415 


2.6225 


-4.89 


0.8959 


-2.34 


1-2 


4.2493 


1.6345 


4.0912 


-3.72 


1.5675 


-2.54 


2-3 


5.6746 


2.3190 


5.4972 


-3.12 


2.2363 


-2.49 


3-4 


7.0961 


3.0126 


7.0275 


-0.96 


2.9845 


-0.69 


4-5 


8.5025 


3.7055 


8.2415 


-3.07 


3.5720 


-2.83 


5-6 


9.7537 


4.3248 


9.6016 


-1.55 


4.2412 


-1.57 



TABLE VI. Magnetic field Bm and at the six 

first degeneracies at e = 0.5 



case, when e > 0.1 already the relation between Bm 
and e is no more linear. On the other hand the spacing 
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Obviously there is a very good agreement, since the 
largest error for Bm is less than 2% and for Em less than 
0.7%. Table VI gives the same but for e = 0.5. Again 
we see a good agreement (error less than 5%). When 
e increases we found that a increases and /3 decreases 
as well as C and our trial wave function becomes less 
accurate, because the double-well nature of the potential 
is less important compared to the kinetic energy. 




2 4 2 4 

B B 

FIG. 7. M values of the ground state depending 
V with increasing e 



Figure |^ describes the situation in the v - B plane for 
m = 0,l,...,10 and different values of e. We notice that 
when V is less than 2 and e is not too large (e < 0.2), 
the situation is similar to the one already discussed, but 
that there is an abrupt change at w = 2 when e is small 
in agreement with the classical analysis. However when 
e > 0.2 the ground state m = is definitely favoured as 
V increases. 




FIG. 8. Comparison for different v of the energy for 
m — 0,1, 2, 3, 4, 5 depending on B with e = 0.1 

Figure ^ shows the energies for the first five m values 
computed with three different v: two-dimensional (v = 
0), spherical potential {v — 2), and v = 3. We can see a 
new crossing between the m — and the other m levels 
when V becomes larger than 2, although this does not 
concern the ground state. 



V. BOUNDS ON THE CRITICAL FIELD IN THE 
TWO DIMENSIONAL CASE 

One might desire to get rigorous upper and lower 
bounds on the critical fields. One possible approach 
would consist in getting upper and lower bounds on the 
ground state energies Em- Whereas we have seen that 
one can obtain very good variational upper bounds, it is 
rather difficult to get good lower ones. In order to test 
these results, we analysed only the two-dimensional case. 

First we want to obtain conditions under which m = 
is the ground state. Using the inequality 



/2 



2P 

X a'- a 

valid for any x and a positive, we deduce that 

2(em)2 



eo[A] > 
On the other hand 



a 



(31) 



(32) 



eo[A]-eo[A-(fi)^]= 



A-(-)^ 



dA'(r2)o(A') 



<(^)'(Oo[A-(^)^] 



(33) 



since (r^)o(A) is decreasing in A. 
But 



A 



o[A]-- < eo[A] 



(34) 



The scaling relation and the fact that eo is increasing in 
e imply that when 4 < — 1 



eo[A] 



< 



(eo[-2] + l) 



(35) 



Taking now a such that ^ > ^ {m> \) vje get combin- 
ing these inequalities that 



Eo < E„ 



Vto > 1 



if we can find t > -f such that 



B^ 



l + -{t^- — ]}S-2e{t- — ]<0 



B 



(36) 



(37) 



where 6^1 + ^eo[-2] + 1 

In the estimate for 6 we can use our best variational 
upper bound. Inequality (|^ will be satisfied if B is less 
than some value Bq, so that in this range m = is the 
ground state. In order to see when m ^ is a ground 
state, we use the following trial wave function ^p{r) for a 
state with angular momentum to'. 



m' > m 



(38) 
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where ipmif) is the exact ground state wave function for 
the state with angular momentum m. An integration by 
parts shows that 



dr r r^'"'-") + 2(m' - m)r2("'-™)-i V^V'm] 
drr2("'-™)^„(rV'™)' (39) 



(40) 



Therefore if we use the fact that 



■{rtl^'mY = [Vm{r) - eyn]tp„ 



We see that 

D POO 

dr r [e V' + 14.' (r) V''] = e™ / dr r 

/•OO 

+ e^2m'{m' -m) dr r^f'"'-™)- V™ (41) 

"'0 

and we conclude that 

/„2(m'-m-l)\ 

em' <em + e'2m'{m' - m) \ (42) 



(r-2(™'-™)). 



In particular 



ei < eo + 2e^ 



1 



(43) 



If we have a lower bound c on (r^)o then we see that 

El < Eo (44) 

if 

2e 



B > 



(45) 



We can use for the lower bound c the one given in equa- 
tion M) 



(46) 



which is satisfactory when B is not too large, but which 
becomes negajtive for large B. We can repair this by 
using the facta that if / is an increasing function of r, 
its expectation value in the ground state is lowered by 
adding to the potential a new increasing potential. We 
can find a useful comparison potential 



W = air^ + a2r^ + asr^ > V 



(47) 



which has a ground state wave function of the form 

^ = e'''r'-b2r' 62 > (48) 

SO that {r'^)w can be computed explicitly for this poten- 



tial and we can take c = {r'^)w in equation (45), which 
gives a more satisfactory result for large B. 



In any case we see that the state m = 1 if favoured over 
the state m = if B is larger than some value, and by 
continuity there must exist a field for which both states 
have equal energy. But in order to prove that the ground 
state is m = 1 when B is in some range requires to show 
that E„i > El Vm > 2. For this purpose let us consider 
m as a continuous parameter. Then 



dE„ 



^ 2e2m(4)m - eB 



^ 2,,. -^^ (49) 

om r^ 

If we can show that > for all to > 1, then we will 
have shown that Em > Ei. When to > 1 we have 

1 



> 



and 



(eto)^ 



+ {r^)l, + \{r^)m < er. 



(50) 



(51) 



In order to get a variational bound on e„i we can use the 

2 

trial wave function ip ~ j.m^-ar ^ ^ff]l[Q]l gives 

1, 



< 



:V, 



m+2{Xm+2) 



m + 2 

where Xm is the value of x which minimises 

(eto)^ 



Vm{x) 



+ x^ + Xx 



Noting that equation ( |51| ) implies that 

„2\ 



{'^ )m — -^m ~t~ \j ^rn ~^ ^rni^^m^ 



(52) 



(53) 



(54) 



one can see by combining equations (^9|), (pO|), ( p2| ) and 
(H) that Em > El for all m > 2 if 



< 



1 



with 



= \\\ Xi + \/Vi+2m{Xi+2m) " Vi{Xi) 



(55) 



(56) 



which implies that B should be less than some value. 
We give in the table VII some numerical values for the 



bounds obtained by these methods. 



e 


Ef) < El 


-Bo-i 


El < Em 


Bl-2 


El < Eo 


0.01 


0.0 - 0.005 


0.011 


0.0 - 0.026 


0.030 


0.022 - 


0.05 


0.0 - 0.024 


0.054 


0.0 - 0.124 


0.163 


0.146 - 


0.1 


0.0 - 0.047 


0.121 


0.0 - 0.221 


0.359 


0.366 - 


0.2 


0.0 - 0.088 


0.340 


0.0 - 0.364 


0.878 


1.788 - 


0.5 


0.0 - 0.189 


1.610 


0.0 - 0.609 


2.745 


2.834 - 


1.0 


0.0 - 0.310 


3.686 


0.0 - 0.826 


5.846 


4.277 - 


2.0 


0.0 - 0.469 


7.816 


0.0 - 1.066 


11.910 


8.141 - 



TABLE VII. Results of the inequalities on the energies and 
values B,„ 
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They show that whereas the range of values of B for 
which Eq < El and Ei < Eq is reasonably estimated 
for e ^ 0.1, there is no range of values of B for which 
our bounds show that m = 1 is the ground state except 
when e is very small (0.01) But in this range WKB works 
perfectly well. Obviously we have too poorly estimated 
the effect of the kinetic energy and that of the centrifugal 
barrier. Numerical computations for example show that 
the replacement of (^)i by j^pji is not appropriate when 
e or B are too large. 

In conclusion, even in two dimensions improved rigor- 
ous bounds on the critical values of the magnetic field 
are needed, and the WKB method for which we have no 
estimate of the error gives the best analytic results. 



VI. CONCLUSION 

It could be of course quite interesting to see an ex- 
perimental verification of these surprising effects of the 
magnetic field. Even though we have found them in the 
case of a double-well, we think that the details of the 
potential do not matter too much. What is needed is a 
potential whose minimum is taken sufficiently far from 
the origin. 

We have thought of two possible fields where one could 
observe such effects. The first one is molecular physics 
where often the dynamics of electrons or protons is mod- 
elled by the motion of a quantum particle in a double- well 
(although admittedly often a one-dimensional one.) If we 
consider the case of the electron in the rotationally sym- 
metric double- well, the smallest value of the critical field 
where the m = 1 and m — states are degenerate is 
about 15 Tesla if we take for the depth of the potential 

1 eV and for the distance to the origin of the minimum 

2 A. For protons the situation is more favourable since a 
field of 5 Tesla can create a degeneracy when the depth 
is kept to 1 eV and the minimum is at a distance of 1.5 
A. Obviously a more detailed investigation is needed if 
one wants to see these unusual effects (like a change from 
diamagnietism to paramagnetism) in molecules. 

The other field is that of Bose-Einstein condensates 
of very cold atoms, which recently has made spectac- 
\ilar progress. If we consider free charged bosons in a 
magnetic field and in a potential ^(7^) one can show 
that there is a Bose-Einstein condensation in the ground 
state in three dimensions, in the limit tq going to infinity, 
for all potentials which have a quadratic dependence of 
r near the origin. Our result supports therefore that free 
charged bosons in their condensate would show a phase 
transition when one varies the magnetic field. This tran- 
sition would manifest itself by jumps of the magnetisa- 
tion at some critical values of the magnetic field. The 
phenomenon would probably persist in a dilute gas of 
charged bosons in a neutralising background. It is how- 
ever probably quite difficult to create such a jellium in 
the laboratory and this remains a challenging task. 
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